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Abstract 



The SUSY-QCD radiative corrections to the T(H + — > tb) partial 
decay width are analyzed within the Minimal Supersymmetric Stan- 
dard Model at the one- loop level, 0(a s ), and in the decoupling limit. 
We present the analytical expressions of these corrections in the large 
SUSY masses limit and study the decoupling behaviour of these co- 
rrections in various limiting cases. We find that if the SUSY mass 
parameters are large and of the same order, the one loop SUSY-QCD 
corrections do not decouple. The non-decoupling contribution is en- 
hanced by tan/3 and therefore large corrections are expected in the 
large tan (3 limit. In contrast, we also find that the SUSY-QCD cor- 
rections decouple if the masses of either the squarks or the gluinos are 
separately taken large. 
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1 Introduction 



Despite its tremendous success in the agreement between the Standard Model 
(SM) predictions and the experimental data, this model contains a variety of 
theoretical problems which cannot be solved without the introduction of some 
new physics. Here we will be concerned about two candidates for physics 
beyond the SM: Supersymmetry (SUSY) and Extended Higgs Sectors. 

The minimal extension of the SM Higgs sector is a two Higgs doublet 
model (2HDM). This model predicts three neutral Higgs bosons, two CP even 
(h°, H°) one CP odd (A°), and two charged bosons (i^ ± ); and includes one 
parameter, tan /3, defined as the ratio of the two Higgs vacuum expectation 
values (vevs) On the other hand, constructing a Minimal Supersymmetric 
extension of the Standard Model (MSSM) implies a two Higgs doublet 
sector, with the same spectrum as in a 2HDM extension of the SM but 
with constrained Yukawa and Higgs self-couplings. To be more specific, the 
MSSM Higgs sector belongs to the so-called type II models, in which one 
Higgs doublet couples to u-like quarks and the other one to d-like quarks 

i- 

Looking for non-SM effects on the Higgs phenomenology of these extended 
models we realize that, while both the neutral Higgs bosons may not be 
easily distinguishable from that of the SM, the discovery of H and the 
determination of its mass and couplings are expected to be a very clear 
signal of physics beyond the SM in the Higgs sector. This is the main reason 
why the search for charged Higgs bosons is one of the major tasks at present 
and future colliders. 

Concerning the present experimental search of charged Higgs bosons the 
situation is as follows. LEP2 has set a model independent lower limit on the 
H + mass, m#+ > 77.4 GeV [Q . At Tevatron, the CDF and DO collaborations 
have searched for H ± bosons in top decays through the process pp — > ti, with 
at least one of the top quarks decaying via t — > H + b, leading to an excess 
of r due to the H + — > t + v decay; they have excluded regions with light H + 
and large values of tan/?, for example if M H + = 100 GeV, tan/3 > 60 is 
excluded [|J. In any case, the allowed region in the (tan/3 — M#±) plane can 
be significantly modified by quantum corrections |J. 

Here we consider the case of a heavy charged Higgs boson, Mh+ > vn t + 
rrib. In this case H + bosons will decay mostly into tb pairs with approximately 
BR > 85% || , depending on the values of tan (3 and M H + . On the other hand 
the two main H + production subprocesses which will provide sizeable cross 
sections at LHC are: gb — > tH + |7j and gg — > tH + b ||. Thus, the importance 
of the H + — > tb decay comes in two directions: it is the dominant decay over 
most of the parameter space and its associated effective vertex H + tb appears 
in the two main production processes. 

Once a charged Higgs Boson is found, a precise determination of its cou- 
plings to SM particles, which are sensitive to radiative corrections, can give 
us indirect information about extra new physics beyond the SM. In particu- 
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lar one may wonder if the H ± is itself an indirect signal of supersymmetry 
and explore if its corrected couplings to SM fermions are like in the MSSM 
or like in a non-SUSY 2HDM. In order to compare the predictions in the 
MSSM and in the 2HDM extension of the SM for these couplings, we assume 
here the most pessimistic scenario where the genuine supersymmetric spec- 
trum is very heavy as compared with the electroweak scale. This situation 
corresponds to the decoupling of supersymmetric particles from the rest of 
the MSSM spectrum, namely, the SM particles and the MSSM Higgs sector. 
In this decoupling limit the SUSY particles can not be produced directly in 
the colliders and we are constrained to indirect searches for SUSY. Here we 
will look for indirect heavy SUSY signals through their effect on radiative 
corrections to the H + — ► tb decay. 

Recent works on the decoupling limit have shown that all the genuine 
heavy SUSY particles and heavy Higgs bosons (H°, A° and if ± ) of the MSSM 
decouple, at one-loop order, from the low-energy electroweak gauge bosons 
physics [|[] . On the other hand, concerning the Higgs physics, the decoupling 
behaviour of the one-loop SUSY-QCD radiative corrections to the T(h° — ► bb) 
partial width have been studied in [|l(|. There it was found that, for fixed 
masses of the extra Higgses (H°, A° and H^), the SUSY-QCD corrections 
from nearly heavy degenerate gluinos and squarks do not decouple in the 
T(h° — > bb) decay width. The SUSY-QCD corrections to the main Higgs 



bosons and top decays have also been analyzed in [|Tl[ . 

In this paper we study the MSSM radiative corrections to the H + — ► tb 
decay width at one loop level and to leading order in a s , and we analyze their 
behaviour in the decoupling limit of heavy SUSY particles. These corrections 
come from the SUSY-QCD [|l|,[13|] and pure QCD sectors |TJ| and are known 
to provide the dominant contributions to the T(H + — > tb) partial width. The 
QCD corrections range from +10% to —50% and the SUSY-QCD corrections 
can be even larger in a large region of the SUSY parameter space. Since a 
2HDM extension of the SM has no SUSY-QCD contributions to H+ -> tb, 
one could use these corrections in order to distinguish between the MSSM 
and the 2HDM. 

Here we analyze the full diagrammatic formulae for the on-shell one-loop 
SUSY-QCD corrections to the H + — > tb decay. We perform expansions in 
inverse powers of the SUSY masses in order to examine the decoupling be- 
havior when these masses are large compared to the electroweak scale. The 
SUSY-QCD corrections depend on a number of different MSSM mass pa- 
rameters, and we will see that the relative sizes of these parameters do affect 
qualitatively the decoupling behaviour. To remain as model-independent as 
possible, we make no assumptions about relations among the MSSM pa- 
rameters that may arise from grand unification or specific SUSY-breaking 
scenarios. We consider the soft-SUSY-breaking parameters and the /x pa- 
rameter as unknown parameters whose magnitudes are all of order 1 TeV. 
We examine in great detail the case of large tan/5, for which the SUSY- 
QCD corrections are enhanced. This enhancement gives rise to a significant 
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one-loop correction, even for moderate and large values of the SUSY masses. 

This paper is organized as follows. In Section ^| we define our notation 
and briefly review the Higgs and squark sectors of the MSSM. In Section [3] 
we review the exact one loop result for the SUSY-QCD corrections to the 
H + — > tb partial decay width. In Section f| we derive analytic expressions for 
these corrections in the limit of large SUSY masses and for several extreme 
squark mixing cases. We also analyze the decoupling of the SUSY-QCD 
corrections for various hierarchies of mass parameters, and compare the ana- 
lytic approximations to the exact one- loop result. Finally, we summarize our 
conclusions in Section ||. 



2 MSSM Higgs and squark sectors 

In the MSSM there are two isospin Higgs doublets containing eight degrees 
of freedom. After the electroweak symmetry-breaking mechanism, three of 
these eight degrees of freedom are absorbed by the Z and W ± gauge bosons, 
leading to the existence of five physical Higgs particles. These consist of two 
CP-even neutral scalar particles h°, H°, one CP-odd neutral pseudoscalar 
particle A° and two charged scalar particles H . Due to supersymmetry, 
the parameters of the Higgs sector are constrained and, at tree-level, it turns 
out that the Higgs masses and mixing angle depend on just two unknown 
parameters. These are commonly chosen to be the mass of the CP-odd 
neutral Higgs boson, M4, and the ratio of the vevs of the two Higgs doublets, 
tan/3 = Vijv\. The Higgs bosons masses at tree level, in terms of these 
parameters, are given by: 

M 2 H± =M\ + Ml 



M 2 - - 

m H°,h° — 2 



M 2 A + M 2 Z ± J (Ml + M 2 ) 2 - AMlM 2 cos 2 2/3 



We choose a convention where the vevs are positive so that < (3 < tt/2. 
The mixing angle in the neutral sector is given at tree-level by: 

tan 2a = tan 20 f. (2) 

M 2 A - M 2 Z w 

In the conventions employed here, —n/2 < a < 0. 

We now discuss the parameters of the third generation squark sector. For 
simplicity, we assume here that there is no intergenerational flavour mixing. 
The tree-level stop and sbottom squared-mass matrices are: 



M ' - { Ml ) ■ < 3 > 
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where: 



M 2 L 


= M| + m 2 + cos 2/3(1/2 


- 2/3 S ^) 


Ml 


Mr 


= M? + m 2 +2/3 cos 2(3 s 2 w M% 




x t 


= A t — /i cot /3 






M ' 2 


= M| + m2- cos 2/3(1/2 


- 1/3^) 


Ml 


Ml 


= M| + - l/3cos2/3 






x b 


= A b — fi tan /? , 






w 


The parameters Mq, 


and My 


are 



(5) 



and sw = sin 6^ 

breaking masses for the third-generation SU(2) squark doublet Q = (t L ,b L ) 
and the singlets D = bji and U = tn, respectively. A b)t are the corresponding 
soft-SUSY-breaking trilinear couplings and /i is the bilinear coupling of the 
two Higgs doublet. The squarks (sbottom and stop) mass eigenstates are 
given by: 



Q2 



where 



COS 6 n 



sin 6 n 



sin 8g cos 9q 



The stop and sbottom mass eigenvalues are given by[]: 



(6) 



(7) 



Mf 

H.2 



Mi 

01,2 



M 2 L + M 2 R ± 



M' 2 + M' 2 ± 



(M- 



M 2 R ) 2 



Am 2 X 2 



{M' 2 



M' 2 Y + Am 2 X 2 



And the mixing angles 9q (q = t, b) are given by: 



cos 26 1 



M 2 



M 2 R 



Mf 



M? 
sin 26 n 



cos 26r 



Ml -M? 

91 92 



M2 2 



Mr 



Ml 

01 



M| 

&2 



(9) 



Concerning the experimental bounds on the SUSY and Higgs masses that 
are relevant for this work, we briefly summarize next the present situation. 
As mentioned before, from combined searches of charged Higgs bosons at 
LEP one gets a general lower limit of M^+ > 77.4 GeV ||, valid at the 95% 
C.L. and for any value of the branching ratio B(H + — > t + v). On the other 
hand, from absence of signals in direct searches at the Tevatron [15|, the 
sbottoms must be heavier than about 140 GeV, assuming that the mass of 
the lightest neutralino Xi is l ess than half the mass of the lighter sbottom. If 



1 Note that in our convention, Mt. 



> M~ t2 and Mr 



> M h 
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m^o = 40 GeV the maximum excluded value for the scalar top mass is about 
110 GeV at the 95% C.L.. For heavier neutralinos, the Tevatron searches 
lose efficiency. In this region the direct searches at LEP [|16|] place a lower 
bound on the sbottom masses of about 91 GeV under the assumption of a 
dominant b — > b\\ decay. For t — Xi mass splittings in the range from 6 to 
40 GeV, i.e., a region not accessible to the Tevatron searches, the lower limit 
on scalar top mass is 83 GeV, independently of the t mixing angle. All these 
lower limit results are at 95% C.L.. Finally, the limits on the gluino mass 
Mg are more model-dependent. If one assumes relations between the gaugino 
masses such that they unify at the GUT scale, then M g is constrained from 
direct searches at the Tevatron to be greater than 173 GeV, independently 
of the squark masses ITi] . 



3 H + — > tb partial decay width 

The H + — > tb partial decay width at tree level is determined by the interac- 
tion Lagrangian describing the H + ib vertex, which is given by [18|: 

C mb = 9^ H + t [m t cot (3P L + m b tan f3P R ]b + h.c. , (10) 
V 2Mw 

where Pl,r = 1/2(1 =F 7s) are the chirality projection operators and V t b is 
the corresponding CKM matrix element. Here we do not consider mixing 
between different families and we set V t b = 1- 

In this work we are interested in the dominant radiative corrections to 
the tree-level partial decay width H + — > tb, which are known to come from 
QCD loops corrections. The leading corrections at one-loop level and to order 
a s in the context of the MSSM were computed in P^-TT4|. They included 
the pure QCD corrections from quarks and gluons and the genuine SUSY- 
QCD corrections from gluinos and squarks. The standard QCD corrections 
were first computed in |14| and can be large, ranging from +10% to —50% 



with respect to the tree level contribution. The SUSY-QCD corrections, 
computed by using a diagrammatic approach in |T^,|13[, have been found to 
be comparable or even larger than the standard QCD contributions in a large 
region of the SUSY parameter space. 

For completeness, and since the SUSY-QCD corrections are our starting 
point in order to study their behaviour in the large SUSY mass limit, we 
have reproduced the results in |L2|,[13]] and we present in the following a short 
summary of the most relevant analytical expressions. 

Following the standard renormalization procedure, there are two kinds of 
contributions to the partial decay width: one coming from the loop diagrams 
and another one coming from the counterterms. Taking this into account, 
we can write the H + — > tb partial decay width at the one-loop level and to 
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order a s in the following wayQ 

T 1 (H + -> tb) = T {H+ -»■ + 2A ioops + 2A CT ) , (11) 

where r\ is the one-loop partial width, T is the tree-level partial width and 
^loops an( j ^ct p rov i(j e corrections coming from the one-loop diagrams 
and from the counterterms, respectively. 

Since the renormalization of the Higgs wave function, Higgs mass, the 
vevs (and hence tan/3) and the parameters g and M w receive no 0(a s ) 
corrections at one-loop, the counterterms contribution originate just from 
the renormalization of the quark wave functions and quark masses. The 
corresponding expression for A CT is: 



D \m t 2 L 2 R J D \m b 2 L 2 R 



(12) 



where: 



D = 


(M 2 H+ 


-m 2 


-m 2 b ) 


(ml cot 2 (3 + ml tan 2 /3) — 4m 2 m 2 


u t = 




-m\ 


-m 2 b ) 


m 2 cot 2 /3 — 2m 2 m 2 , 


u b = 


(M 2 + 


~ ™t 


-m 2 b ) 


m 2 tan 2 /3 — 2m 2 m 2 . 



(13) 

We use the on-shell renormalization scheme as defined in [Tj|. The contri- 
butions from the counterterms can be written in this scheme in terms of the 
bottom and top self-energies as follows: 

+ Uz^ + Uz™ = 5#»W W ) + ^ViU) - *#V 2 ' 

Ht,b) 



2 



££(m 2 ) + 4K) + 2E£(m 2 ) 



m 2 . 



K 2 ) + ^(m 2 )+2S^(m 2 

(14) 



As we have already mentioned, in the present work we focus on the con- 
tributions to A' oops and A CT coming from the SUSY-QCD sector: 

&SQCD = ^SQCD + ^SQCD 

where we use the short notation SQCD for SUSY-QCD. 

The SUSY-QCD contributions to the H + — > tb decay width come from 
diagrams involving the exchange of virtual gluinos (g), sbottoms (b) and 
stops (t). The diagrams that contribute to the vertex corrections and quark 

2 The factor 2 used here is a convention such that the corresponding correction to the 
effective coupling is: 

W P = 9%? b & + A' oops + A CT ) • 
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Figure 1: Diagrams corresponding to the vertex correction and quark self- 
energies that contribute to A S qcd- 



self energies are shown in Fig. [I]. The contribution from the one-loop vertex 
diagrams are given by: 



±s35D = %H t + %H ll . (16) 



where: 



H t = — ^ — ^ - [m t R®Ria* (C u - C 12 ) + m b R^R^*Ci 2 



+M~ g RVR?>C ](ml Ml +1 ml M|, M|, MfJ , 



+M- g R^R { S*C ](ml M 2 H+ ,ml Mf , M|, MfJ , (17) 

R^ are given in eq. (0) and we have used G ab to parametrize the if + b a t b 
coupling. The expression for G ab is given in the Appendix A. 

Finally, the contributions from the SUSY-QCD sector to the quark self- 
energies are: 

mv 2 ) = - 2 -^\r^\ 2 b 1 (p 2 ,M|,My, 



37T 

S|(P 2 ) = Re( 4? 4f )5o(p 2 , M|, 



^W) = -^m 2 Bi{p\MiMi: 



<4 



We have decomposed the quark self-energy according to 

£^) = £{(p 2 ) Z>P L + 4(p 2 ) i'PR + m^^p 2 ), (19) 
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and used the notation S'(p) = <9£(p)/<9p 2 . 

Our notation for the one-loop integrals that appear in the above formulae, 
B , B' Q , Bi, B[, Cq, Cn and Cu, is defined in the Appendix B. We have 



checked that these results are in agreement with the calculation of [12, 13 



4 Large SUSY mass limit in the MSSM 

Here we study the effect of heavy SUSY particles in the SUSY-QCD correc- 
tion to H + — > tb. We derive approximate expansions for this correction in the 
limit of large SUSY mass parameters as compared to the electroweak scale, 
M EW , which corresponds to a SUSY spectrum heavier than the Standard 
Model one. 

To define our expansion parameters, we consider all the soft-SUSY-brea- 
king mass parameters and the n parameter to be of the same order (co- 
llectively denoted by Msusy) and much heavier than the electroweak scale. 
That is, 

M SUSY ~ Mq ~ M v ~ M 3 ~ M~ g ~ p ~ A t ~ A b > M EW , (20) 

where all these parameters are defined in Section |2|. Notice that, with 
Msusy ~ 0(lTeV), this choice will lead to a plausible situation where all 
the SUSY particles in the SUSY-QCD sector are much heavier than their 
SM partners.^ 

We compute the expansions of the SUSY-QCD correction to the partial 
decay width T(H + — > tb) in inverse powers of the SUSY mass parameters, up 
to order M EW /Mg USY , considering that all M H +, M z , M w , m t and m h are of 
order M EW . In this section we present the leading terms of these expansions 
while the complete expressions are shown in the Appendix C. 

We consider here two possible extreme configurations of the squarks mass- 
squared matrices corresponding to maximal and near-zero mixing respec- 
tively. This situation leads us to four extreme cases that will be considered 
in the present work: 

• Maximal mixing in the sbottom and stop sectors (0g ~ ±45°, Of ~ 
±45°) , 

• Near zero mixing in the sbottom and stop sectors (0g ~ 0°, 0j ~ 0°) , 

• Maximal mixing in the sbottom sector and near zero mixing in the stop 
sector (6~ b ~ ±45°, % ~ 0°) , 

• Near zero mixing in the sbottom sector and maximal mixing in the stop 
sector (6- b ~ 0°, 9~ t ~ ±45°) . 

3 Notice that our choice of large [i is not a necessary condition to get all the SUSY-QCD 
particles heavy but it is needed in order to get all the gauginos heavy in the electroweak 
sector. 
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Maximal mixing (9 q ~ ±45° ) arises when the splitting between the dia- 
gonal elements of the mass-squared matrix is small compared to the off- 
diagonal elements, that is \M\ — M R \ -C m 4 |X t | in the stop sector or \M L 2 — 
M R \ <C m |X | in the sbottom sector. This situation leads to small mass 
splitting between the two squarks compared with the masses themselves, 
namely, \M? — M~ 2 | <C |M? + M| 2 | . In order to define properly the large 
SUSY mass expansions in powers of a small expansion parameter, we notice 
that all mass scales should be referred either to Mew or to Msusy- Thus, the 
counting in this maximal case goes as follows: m q X q is of order MewMsusy, 
while the splitting between the diagonal mass-squared terms is of order M EW , 
so that the mass-squared eigenvalues can be written as: 



M~ 

51,2 



Mf ± A 2 



or /and 



M~ 

01,2 



Mj±A 2 



(21) 



where we have defined: 



M 2 s = U M i+ M l 



Ml 



UK 



+ M?' 

b 2 < 



(22) 



A 2 



m t \X t 



m b \X b 



1 + 



1 + 



(M- 



Ml? 



im 2 t Xf 



(M' L 2 



8m 2 X 2 



M'if 



(23) 



and M L , M R , M L and M R are defined in eqs. (|5]). Here M| and M| are 
of order Ml USY while A 2 and A 2 are of order MewMsusy ■ The quantities 



Ml- 



and 



'2 



are therefore small and of order Mew/M, 



mtx t m b x b — — — ^ ^ ^ ^mw/^susY- Expan- 

ding the expressions for the mixing angle in terms of these small parameters, 
we obtain for 9 q ~ 45°: 



cos 29 1 ~ 



Ml 



Ml 



2m t X t 
sin 29, 



sin 26^ ~ 



cos 2^ ~ 



M'l 



(Mj 



Mr 



2m b X b 
Ml? 



8mjXf 

{Ml - M'j? 
%m\Xl 



(24) 



where <Jx q = sgn(X q ) with q = t,b. 

Near-zero mixing (9 q ~ 0°) arises when the splitting between the diagonal 
elements of the mass-squared matrix is large compared to the off-diagonal 
elements, that is |M| — M R \ 3> m t \X t | in the stop sector or \M'? — M R \ ^> 
m b \X b \ in the sbottom sector. This is the case usually considered in the 
literature, because M L (M L ) and M R (M R ) depend on two different (and 
a priori independent) soft-SUSY-breaking parameters Mq and M v (Mf)), 
respectively. The counting in terms of Mew and Msusy is such that (M| — 
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Mr) (or [M'l — Mr)) is of order M SUSY while m g X g is of order M EH /Ms;ysy. 
The mass splitting between the two squarks will be of the same order as the 
masses themselves, that is |M? - M|| ~ 0(|Af| + M?|) ~ C^M^y). As 
in the maximal mixing case, all mass scales should be referred to Mew or to 
Msusy in order to define the expansion parameters. In this case we write 
our results in terms of the physical stop and sbottom masses and the mixing 
angles which are given by: 

A ^=w^m' < 25) 

91 12 

2m 2 n X 2 n 

cos2 ^ 1 -(^%' (26) 

4.1 The decoupling regime in H + — > tb decay 

Now we study, both analytically and numerically, the SUSY-QCD correction 
to T(H + — > tb) in the large SUSY mass parameters limit. We consider here 
all the possible extreme configurations of the squarks mass-squared matrix, 
namely the four cases mentioned above. We will study also the individual 
decoupling of the various SUSY particles, that is, decoupling of gluinos and 
squarks separately. In this analysis we use the leading term of our expansions 
for A S qcd while the complete expressions, up to order M EW /M SUSY , are 
given in the Appendix C. 

First we consider maximal mixing in both the stop and sbottom sec- 
tors. In this case, and in the large SUSY mass parameters limit defined in 
eq. (pOD , we obtain M| ~ Mr and M' 2 ~ Mr, up to corrections of order 
M EW /Mg USY . Since Ml and M L are determined by the same soft-SUSY 
breaking term Mq, we also have M| ~ M 2 . From these relations and from 
eq. fl23|) we obtain M| ~ Mf , up to corrections of order M EW /M SUSY . 

Taking this into account we expand the SUSY-QCD contributions to the 
H + — > tb partial width and we find: 

^sqcb = -^^|Vn/? + cot + 0(M 2 EW /M 2 SUSY ) (27) 

where R = Mg/Ms- The expression for fi(R), which arises from the expan- 
sions of the loop integrals, is given in Appendix B and it is normalized so 
that /i(l) = 1. The complete formula for A S qcd, including the Oi Ai2 EW ) 

M SUSY 

terms, is given in Appendix C. 

We can see from eq. (|2T| ) that, taking all SUSY mass parameters arbitrari- 
ly large and of the same order, Asqcd leads to a non-zero value. That is, the 
SUSY-QCD correction do not decouple in the large SUSY mass parameters 
regime. This can be seen clearly, for instance, in the simplest case of equal 
mass scales, /i = M g = Ms, where fi(R) = 1, leading to 

&sqcd = -^(tan/3 + cot (3) + 0(M 2 EW /M 2 SUSY ) } (28) 
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Figure 2: Non- decoupling behavior of A S qcd with Mq = M v = = M g = 
A b = A t = n = M S usy an d for different values of tan p. Exact one-loop 
result (solid lines) and the expansion given in eq. ft2l[) (dashed lines) are 
plotted for comparison. 



which shows that the first term in the large Msusy expansion is indeed of 
O{M% w /M SUSY ). This leading contribut ion, when considered in the large 
tan P regime and expressed in terms of an effective coupling of H + to bi 



is in agreement with the previous results of refs. |20| that were obtained in 
the zero external momentum approximation by using an effective Lagrangian 
approach. 

Notice also that this non-decoupling term is enhanced by a tan f3 factor 
and therefore a numerically large SUSY-QCD correction in the large tan/3 
limit is expected. This non-decoupling behavior is shown numerically in 
Fig. ^. In this figure we plot the exact one-loop result (solid lines) and the 
leading term of the large SUSY mass expansion in eq. (p7|) (dashed lines) for 
Mq = Mfj = Mf) = M g — A\ — A t — fM — Msusy and for different values 
of tan p. For all the numerical analyses, in this work we fix Mjj+ = 250 
GeV. We can see in this figure that for large SUSY mass parameters (say 
Msusy > 300 GeV), A S qcd tends to a non-zero value whose size is larger for 
larger tan P values. Indeed the correction can be quite large for large tan P 
even for quite heavy SUSY particles. For example for Msusy = 1 TeV and 
tan/3 = 30 we get Asqcd — —35%, and it grows linearly with tan/3. 

From the numerical comparison in Fig. |^ between the exact and our 
approximate result of eq. (|27|), we can conclude that this expansion, with 
just the leading term, is a good approximation for large enough SUSY mass 
parameters. It is clear also that as tan/3 grows the agreement between the 
exact and approximate results becomes worse at low Msusy- However, we 
can conclude that the agreement is still quite good even for tan P as large as 
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Figure 3: Behaviour of the contributions to Asqcd of 0(M% w /M^ USY ) for 
Mq = Mfj = M D = Mg = A b = A t = n = M S usy o,nd for different values of 
tan p. 



50 whenever M SU sy > 300 GeV. 

For lower values of Msusy the next to leading corrections, i.e. the co- 
rrections of order Mg W /Mg USY which are given in Appendix C, begin to be 
important. We have estimated the size of these corrections for several choices 
of the parameter space and we have found that they are of order of a few 
percent. Besides, they grow with the \x parameter and with tan/? and are 
almost independent of the value for the trilinear couplings, A t> b- We show in 
Fig. |] the numerical results for the particular case Mq = M v = Mfj = M g = 
Ab = A t = fi = Msusy and for several values of tan/3. We can see clearly 
that their size is always small as compared to the leading non-decoupling 
term. For other choices of the parameters similar results are found. 

In eq. (p7|) we can see that, in the approximation of large SUSY mass 
parameters, changing the sign of Mgfi simply flips the sign of the SUSY- 
QCD correction. In most of this work, and for the numerical analysis, a 
positive sign for M g [i has been chosen. 

In Fig. f| we show the exact SUSY-QCD correction as a function of \i 
(left) and as a function of tan/? (right). In the left panel we plot Asqcd for 
differents values of tan/3 and fixing Mq = M v = M^ = M g = Ab = A t = 1 
TeV. We see that the size of the correction grows linearly with /i and notice 
that by changing the sign of \x the sign of the correction changes. In the right 
panel we plot the SUSY-QCD correction as a function of tan f3 for the choice 
Mq = M v = Mjj = M g = A b = A t = [i = Msusy and for different values 
of Msusy- We see explicitely that for large Msusy the size of the correction 
grows linearly with tan (5. 

In comparing the various contributions to Asqcd we find that the domi- 
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Figure 5: Comparison of contribution coming from bm^jm^ (lower dashed 
line), SZb (upper dashed line), 5m t /m t ( lower dotted line), SZ t (upper dotted 
line) and vertex contributions (long-dashed line) for different values of tan (3. 
Total correction is also plotted (solid line). 



nant contributions come from the renormalization of the bottom-quark mass 
and wavefunction. In Fig. |5] the total SUSY-QCD correction (Total), the 
contributions from the bottom-quark mass (Smb/rrib) and wavefunction (SZb) 
counterterms, the contributions from the top-quark mass (Smt/mt) and wave- 
function (SZt) counterterms and the vertex contributions (Vertex) are shown 
separately. We can see that the contributions from Smb/mb + SZb dominate 
Asqcd at large tan/3. 

For the three cases left, near-zero stop and sbottom mixing, near-zero 
stop-maximal sbottom mixing and near-zero sbottom-maximal stop mixing, 
we have found a leading term in Asqcd that is very similar to the one in 
eq. Q27p. The general expression is: 

Asqcd = -^-|^(tan/3 + cot (3)F mix + 0(M 2 EW /M 2 SUSY ) , (29) 
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where M| and F m i X depend on the case: 

• Near-zero stop and sbottom mixing: 

M 2 S = M? and F mix = h(R u R 2 ) where Ri = M~ g /M^ and 



h 

M fi ~M; ~il%,i=l,2 



Near-zero stop-maximal sbottom mixing: 

M 2 S = M 2 S and F mix = f/^, i^) + ^f\(R b ) where R t2 = M g /M h 
and R b = M~ g /M s . 



Near-zero sbottom-maximal stop mixing: 



M 2 S = M 2 S and F mix = %A(R t ) + %h(R t ,R b2 ) where R b2 = M g /Mr 



and R t = M- g /M s . 

The functions fi(R) and fi(Ri,R 2 ) are defined in the Appendix B and are 
normalized such that /i(l, 1) = /i(l) = 1. M s and Mg are defined in 

eq. (E2l). The corrections in eq. (E9|) of (9( , ) are given explicitely in 

1 V 1 v M SUSY 

Appendix C. 

In eq. ( P^ ) we see that the non-decoupling contribution to A S q C d appears 
in all the cases so that we can conclude that the SUSY-QCD corrections do 
not decouple in the large SUSY mass parameters regime. We also see that 
Asqcd grows again linearly with tan/3, thus we expect large corrections in 
the large tan j3 limit for all the cases. 

Now we study the decoupling behaviour of the SUSY-QCD corrections as 
individual SUSY particles become heavier than a reference SUSY mass scale 
so that there is a hierarchy among the various SUSY particle masses. We 
will consider two cases: large gluino mass with maximal mixing in stop an 
sbottom sectors and large squark masses (stop and sbottom) with maximal 
mixing in both sectors. 

First, we examine the case of a very heavy gluino compared with the rest 
of the SUSY spectrum and with maximal mixing in sbottom and stop sectors. 
In this case M g 3> Msusy 3> Mew- We evaluate the SUSY-QCD correction 
when the gluino mass is taken very large while the rest of the SUSY mass 
parameters remain fixed to a reference value of the order of Msusy (larger 
than the electroweak scale). Our result for A S q C d is: 



isQCD = + cotP) (l _ log I) ^ + 0( _J^_) . (30) 



In eq. (|30|) we see that the SUSY-QCD correction decouples when we take 
the large gluino mass limit while the rest of the SUSY mass parameters (/i 
and Mg) remain fixed at some scale of the order Msusy- We also notice that 
this decoupling is very slow: A S qcd falls off when increasing M g as 9 . 
We see again that the correction is enhanced by a tan (3 factor, so we expect 
large SUSY-QCD corrections in the large tan/3 limit. In Fig. E we plot the 
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Figure 6: Behaviour of Asqcd in the large M g limit with fixed Mq = My = 
M D — Af, — A t — n — 1 TeV and for different values of tan (3. 



exact result for Asqcd as a function of M g and keeping the other SUSY mass 
parameters fixed at 1 TeV. We can see the slow decoupling with the gluino 
mass and the numerically large correction even for very heavy gluino; for 
example \iM g = 1 TeV and tan (3 = 30 we find Asqcd ^ -40%. Notice that 
the size can be so large that the validity of the perturbative expansion can be 
questionable. We refer the reader to refs. [^(J where this subject is studied 
in more detail and some techniques of resummation for a better convergence 
of the series are proposed. 

Next we analyze the large squark (stop and sbottom) mass limit with 
maximal mixing in both sectors. We consider the case with Ms = Ms ^> 
M g — A t — Ab — [A Mew- If we expand the eq. ( |27|) in the limit of large 
squark masses we obtain: 



limit if the rest of the SUSY mass parameters (M g and fi) remain fixed at 
some scale larger than the electroweak scale. This decoupling is faster than 
the previous case since Asqcd is proportional to the inverse of M|, that is 
the scale of the squark masses. We also see that the correction is enhanced 
by a tan (3 factor. Finally we have also confirmed this decoupling behaviour 
of squarks numerically, as can be seen in Fig. [7], where the explicit result of 
Asqcd, as a function of Msijsy = Mq = M D = Mjy, is shown for several 
values of tan (3. 
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Figure 7: Behaviour of Asqcd in the large squark masses limit with fixed 
Mg = A b = A t = fi = 1 TeV and for different values of tan (3. 

5 Conclusions 

In this work we have looked for indirect SUSY signals through the effect 
of radiative corrections to the H + — > tb decay in the limit of large SUSY 
masses. We have focused on the well known large SUSY-QCD corrections to 
the T(H + — > tb) partial decay width. We have studied in detail the one-loop 
SUSY-QCD corrections to T(H + — > tb) when the SUSY particles are heavier 
than the SM ones. In order to understand analytically the behaviour of these 
corrections in the large SUSY mass limit, we have performed expansions of 
the one-loop partial width that are valid for large values of the SUSY mass 
parameters compared to the electroweak scale and for various choices of the 
mixing in the squarks sector. 

We have shown that for large SUSY mass parameters and all of the same 
order, the SUSY-QCD radiative corrections do not decouple in the H + — > tb 
decay. In other words, if there are heavy squarks and gluinos, with masses 
of the same order and much larger than the SM particles, they will produce 
a non-vanishing contribution, via SUSY-QCD radiative corrections, to the 
potentially measurable partial width T(H + —> tb). We would like to remind 
that the conclusion about decoupling or non-decoupling of SUSY particles 
in the MSSM is not inmediate, since the decoupling theorem does not apply 
to this case, because the MSSM is a theory involving chiral fermions and the 
spontaneous breaking of the electroweak symmetry. Therefore, the SUSY 
decoupling behaviour must be explored case by case. In this direction it has 
been shown by an explicit calculation |J that there is decoupling of heavy 
SUSY particles in observables with external electroweak gauge bosons. 

We interpret the non-decoupling effect found in the present work (as well 
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as the one found in |I0| ) as follows. If we start from the full MSSM as the 
theory valid at large energies and integrate out the heavy SUSY spectrum, 
we are left with a low energy theory, valid at the electroweak scale, which 
contains the Standard Model particles plus two full Higgs doublets. The 
interactions of this effective low energy theory are a priori unknown, but 
they can be derived by explicit integration of the heavy SUSY particles [^] . 
The Higgs-fermion-fermion interactions of the MSSM are like those of a two 
Higgs doublet model (2HDM) of type II. In a 2HDM of type II there are 
some restrictions on the allowed couplings, which in the MSSM case are a 
consequence of supersymmetry. Once the heavy SUSY particles are inte- 
grated out, new low energy effective Higgs-fermion-fermion interactions with 
no restrictions anymore can be generated since SUSY is a (softly) broken 
symmetry. Indeed, one expects to find the most general 2HDM of type III, 
where both Higgs doublets can couple to both u and d type fermions. How- 
ever, to show this one must perform explicitly the computation and find out 
the specific values of the generated effective couplings [02 . 



We also emphasize here that these non-decoupling SUSY-QCD correc- 
tions could give us an indirect signal of supersymmetry at present and future 
colliders even for a very heavy SUSY spectrum at the TeV scale. In par- 
ticular they can provide some clues in the indirect search of a heavy SUSY 
spectrum in the LHC [|23| and next linear colliders. 

We have also examined in this work some special cases in which there is a 
hierarchy among the SUSY mass parameters. In the case of maximal squark 
mixing with Ms = M$ large and the other SUSY mass parameters of order 
a common mass scale M (chosen such that Mew <M< Ms), the SUSY- 
QCD corrections decouple like M 2 /M|. In addition, we have examined the 
case of a large gluino mass with the other SUSY mass parameters of order 
a common mass scale M (chosen such that Mew <C M -C Mg). In this 
case we have found that the SUSY-QCD corrections decouple more slowly, 
like (M/Mg) log(M|/M|). This indicates that sizeable indirect signals from 
gluinos as heavy as several TeV can be obtained in the H + — > tb decay. 

Finally we have seen that, in all cases, the corrections are enhanced by 
a tan (3 factor and therefore large SUSY-QCD contributions to the H + — > tb 
decay are expected for large values of tan /3, even for a very heavy SUSY 
spectrum. 
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Appendices 



A H t b interaction Lagrangian 

The H + 1* b and H~ t b* interaction terms are given by the following interac- 
tion Lagrangian [flj: 



£ H ± tl = — f^ M H \9ll b* L t L + 9rr b* R t R + g LR b* R t L + g RL b* L t R ^j + h.c. , 

(a.i; 



V2M W 

where: 



gi,L — sm 2/5 — (mj cot (3 + m\ tan 0) , 
gRR = -m t m b (t&n(3 + cot/?) , 
9lr = -m b {n + A b tanP) , 

9rl = —m t (fM + A t cot j3) . (A.2) 

In the mass eigenstate basis for the squarks, these interaction terms can be 
parametrized as follows: 

where 

Gab = Rla*Rlb9LL + R2I* R^b 9RR + R^Rl^LR + Rl a *^2b9RL ■ (A.4) 

B Large mass expansion of loop integrals 

Here we define the notation for the two- and three-point integrals that appear 
in eqs. (O) and ([18]) and give formulae for their expansions in inverse powers 
of the SUSY mass parameters. 



We follow the definitions and conventions of [^TJ. The integrals are per- 
formed in D = 4 — e dimensions and the divergent contributions are regular- 
ized by A = - — 7s + log (47r), with the corresponding inclusion of the energy 
scale given here by yUo- 

The two-point integrals are given by: 

4— D f d k -flj i r n , , o o 2\ 

"° J {2,r r ~ ml] [{k + qf - ml] = 16^ {B ° ; Q Bi} iq ; mi ' m2) • 

(B.5) 

The derivatives of the two-point functions are defined as follows: 



d 



(B.6) 
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Finally, the three-point integrals are given by: 



d D k {l;k"} 



(2tt) d [A; 2 - m\\ [(k + Pl ) 2 - m 2 } [{k + Pl + p 2 ) 2 - m 2 ] 

{C ;p1C n + p%C 12 } {p\,pl ) p 2 ]ml ) ml ) ml) , (B.7) 



16vr 2 

where p = —p\ — p 2 . 

We next give the large Msusy expansions of the one-loop integrals in the 
four different cases that we have mentioned in section [|. 

B.l Maximal sbottom and stop mixing 

The results of the large Msusy expansions of the loop integrals are as follows: 

C (m 2 ,M 2 H+ ,m 2 b ;M 2 ,MlM 2 b] ) 

~ ~2Ml fl{RuRb) ~ ^Ml h{RuRb) ~ ^Ml h{RuRb) ~ i^Ml h{RuRb) 



A 4 A 4 • -A 2 A 2 

^f 7 (R t ,R b ) - *f 8 (R uRb ) - (_i)^_J_| /9(i2t)i2fc) 



CuimlMlt^mlMlMlMf.) 

' f 10 (R t ,R b ) + (-ly^xMRuRb) + (-l) j ^rfi2{Rt,R b ) 



~ 3M| JiUV ™ ° ! v ' 8Ml J11K ™ °' K ' 8M 4 
C 12 (m 2 t ,M 2 H+ ,mlM?MlM? 



1 . . ^ A? . , A? _ 



f 13 (Rt,R b ) + {-l) l -rh-JiA{RuR b ) + (-lYznrfaMBt,* 



~ QM 2J16y r ' °' y ' 2AMl J1 ^ r ' °' y ' 12M 4 

M? A 2 
BoK; Ml Ml) ~ A - Zn-^ + ^i(^) + (-1)^/1^) 

m 2 A 4 . m 2 A 2 

+ 6iV% /2 W + 6M| /2 W + (_1)4 12Mf (2/s ~ U){Rq) ' 

1 M~ A 2 

B l (m 2 ; Ml M 2 ) ~ --(A - Jn^JL) + £ 2 (i? 9 ) - (-l)'^/^) 

hW-JkrAVz-hW,)- (B.8) 



12M 2 ' v 97 24M 4 



Where Mj and A 2 are defined in eq. (^) and the functions /j and gi are 
given in terms of the ratios R t = M g /Ms and R b = Mg/Ms as follows: 



fi(Ri, R 2 ) 



21 lln^ + RllnRl-RllnRj 



(1-R 2 )(1-R 2 )(R 2 -R 2 ) 



R 
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f 2 (Ri, R2) 



12R 2 



[R\ - Rt 



R\ — 3R 2 R 2 
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R 2 

+ R\ln-± - R\R\lnR\ - 3R\R\lnR\ - RfR 2 2 lnR 2 2 - ARfR^lnRj 



+ 2R 4 2 lnR\ + 2RlR 4 2 lnR\ + 2R 4 1 R 4 2 lnR\ + 3R 4 R 2 2 lnR 2 2 
- R & 2 lnR\ - R\RllnR 2 ] 



f 3 (Ri,R 2 ) — fi(R\ ) R2) 



\2R\ 



[3R\ - 2R\ 



(l-R 2 2 ) 2 (l-R 2 ) 2 (R 2 -R 2 y 
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d2 d2 d2 d2 

^ t-)2; • fl 2 1 nn2n2i - fL 2 , n n4D4/ - ft 2 , OD 4) -"^ 



- i2?Zn-| - i^Zn-f + 2i&R?Zn-f + 2R%R\ln^ + 2i?|/n , 



fl5(Rl, R2) 



QRl 



[ — r\ + i?^ — r 2 — ?>r\r\ 



Rl{l-R*Y{l-R{){Rl-Riy 

+ R\R 2 + 3i?^i?2 "I" 4i?2 — 4i?^i?2 — 3i?2 "I" 3R 2 R 2 — 2R 2 R®lnR 2 

R 2 

+ QR 4 R 4 lnR 2 2 - QR 2 R\lnR\ + 2R 8 2 lnR 2 2 + 2R 2 2 R\ln-^ 



QR A 2 R\ln^ - 2R%R\lrM 
R^ Ri 



f 1 (R)=f 1 (R,R) 
3 



[1 _ R 2 + #2^2] 



/ 2 (i?) 

/ 4 (i?) 
<7l(#) 
92(R) 



(1-R 2 ) 3 
2 

{l-R 2 Y 
2 

(1 - R 2 Y 
1 

1 



4(1 - i? 2 ) 2 



(1-R 2 ) 2 
[l - R 4 + 2R 2 lnR 2 ] 

[2 + R 6 - QR 4 + 3i? 2 + QR 2 lnR 2 ] 

[1 - 6R 2 + 3R 4 + 2R 6 - mHnR 2 ] 

[l-R 2 + R 2 lnR 2 ] 

[-3 + AR 2 -R 4 - AR 2 lnR 2 + 2R 4 lnR 2 ] 



(B.9) 



Note that in the special case M„ 



Mc 



Ms, R t = R h 



1, the functions 



above are normalized as 1) = 1, /i(l) = 1 and = 0. Notice also 
that we use here the same notation as in ref [10 . 



B.2 Near-zero sbottom and stop mixing 

The loop integrals are expanded as follows: 

1 



C (m 2 t ,M 2 H+ ,mi;M 2 ,Ml,Ml 



2M? 



fi (Ru , Rb 



2AM 4 ^ Ru ' Rb ^ ~ 2AMT^ Ru,Rbj " > ~~ 2AM?^ Rt ' ,R 
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C xl K 2 , , m 2 ; M| , M| , M?. ) ~ / 10 , ) , 

Cia (m 2 , M 2 H+ , m 2 b ; M| , M| , JWg ) ~ ^ /i 3 , Jfc, ) , 

M 2 m 2 
BoK; M 2 , M|) ~ A - In-f + 9l (R qi ) + jr&hiR*) , 



1 M? _ /;/: 



i 2 

^(m 2 ; M 2 , M|) ~ -±(A - Inj^) + g 2 (R qi ) - ^/s^) • 



(B.IO) 

Where R qi = M g /M qi (i — 1, 2) and the functions /, and ^ are the same as 
in Section B.l. 

B.3 Maximal sbottom mixing and near-zero stop mixing 

The loop integrals are expanded as follows: 

C (m 2 ,M^ + ,m 2 ;M 2 ,M 2 ,M 2 ) ~ -JL/^, R b ) 
m 2 „ , „ „ s M 2 , „ , „ „ s m 2 



* jf2{Rti,Rb) - OA f^f3i R t v R b) - OA r74 h(Rti,Rb) 



2AM? J y ' 24M? Joy ' 24M 

Ej Lj_ L<i 

A 2 A 4 
-(-^y^^MRkjRb) - 12 ^_ 6 /8(-R* l ,-Rb) , 

C n (m 2 ,M| + ,m 2 ;M 2 ,M|,M|) ~ /^.ifc) + (-1)^^^, 

C 12 (m 2 t ,M 2 H+ ,m 2 b ;MlMl,Ml) ~ J^/ 13 (i* ti , + (_1)^ ^ /l5 (i^, 

(B.H) 

Here = M g /M~ u (i = 1,2); R b = M g /M s and the functions ft are the 
same as in Section B.l. For the two-point integrals Bo } i(ml, M~, M~) and 
B 0) i(m 2 , M~, M|) we use the expressions given in (|B.8|) and flB.10| ) respec- 
tively. 

B.4 Maximal stop mixing and near-zero sbottom mixing 

The loop integrals are expanded as follows: 
C (m 2 ,M 2 H+ ,m 2 b ;M 2 ,Ml } M 2 b ) 

~ ~2Ml fl{RuRb > ] ~ i^tl f2{RuRb ^ ~ 2^Ml h{RuRb ^ ~ MII hiRuRb ^ 
A 2 A 4 
~ { ~ 1)l 6Ml f5{RtjRb ^ ~ 12Ml fr{Rt,Rbj) ' 

24 



C n (m^M| + X;Mf,M|,M|) ~ -L /l0 (i^ + (-iy^Lf n (R t ,R bj 
C 12 (mlM 2 H+ ,m 2 b] MlMl,Ml.) ~ ^f 13 (R t ,R bj ) + (-iy-M_f 14 (R u R b 



(B.12) 



where R b . = Mg/M b . (i = 1, 2); R t = Mg/Ms and again the functions /j are 
as in Section B.l . In this case for the two-point integrals Bo^im 2 , M~, M~ ) 
and B 0j i(m 2 , M~ , M|) we use the expressions given in ( B.10|) and ( B.8| ) re- 
spectively. 



C Complete expressions for Asqcd to order 
Mew/Msusy 

Here we give the complete expressions for the expansion of the SUSY-QCD 
correction to T(B+ -> tb) up to order M 2 EW j M 2 SU SY . From eqs. (H), (||) 
and (|1^) the SUSY-QCD correction can be rewritten as: 



&SQCD = jj(K t + H t ) + j}(K b + H b ) 



(C.13) 



where we have defined: 



K t = ^ + l -5Zl+ l -5Z^ 
m t 2 L 2 R) 

K b = ^ + Uz\ + Uz\. (C.14) 
m b 2 2 

In the following we give the expressions for K t , B t , K b and B b up to 
order M% W /M^ USY in the four different cases that we have considered in 
this paper. 

C.l Maximal stop and sbottom mixing 

The results for K t , B t , K b and B b up to order Mg W /Mg USY are the following: 



3tt 



2 



M 5 ( ^ /i(i? t ) - ^ (2/ 3 - / 4 )(i4)J - - lnj± + 9l (R t ) - 9l (R 



+ g 2 (R t ) - g 2 (R b ) - ^ (2/ 2 - / 3 )(^) - ^ (2/ 2 - f 3 )(B b 
m 2 X 2 m 2 X 2 

+ 2Im| (4/2 - 2/3 + /4)(jR<) - s| (4/2 - 2/3 + /4) ™ 

M 2 - M 2 M' 2 - M' 2 
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H* = 



2a f 



M^f^R,) + ! |gr( 2 /io - fi3)(Rt,R b ) 



3n nil cot 8 L ^ ^ 

+ ~Fr7r" m tXtf^ > {Rti Rb) + affi 1 hs(Rt, R-b) + /.^ "^^(^t, Rb) 

m t m b X b g LR m t m b X t g LR 
2AMI Wi2-2/i 5 )(^,^) 24M| /l4 ^'^ 

M- g m t m b X t X b g LR m 2 t X t g RL 

12Mg M " &) + 24M| (/l4 ~ 3/l1 ^' ^ 



M g m 2 X 2 g RL M g m 2 X 2 g RL 

V2Ml fl[Ru Rb) 12M| h[Ru Rb 

M g (M' 2 - M'j)g RL M g (M 2 -M 2 R )g RL 

+ ^ HRuRb) ^ HRuRb ^ 

M ~ 9 (f| fl{Rb) " (2/s " /4)(jRt) ) + ^ n f " 9l{Rt) + 9l{Rb) 

1 2 

m in 

- ff2 (i? t ) + g 2 {R b ) - ^ (2/ 2 - h){R b ) - ^ (2/ 2 - / 3 )(it! t ) 



gg(4/ 2 - 2 /3 + /«)(*) + gg(4/ 2 - 2/ 3 + *)(*) 



mlX? 



+ 



Mr — M 



24M 4 . 

M' r 2 - m: 2 



H h = 



2a, 



1 



^«/ 1 (i2 t , J R 6 ) + ^^(2/io - fis)(Rt,Rb) 



3n mi, l an .'^ L 

+ ^Mf mtXtm - Rb) + w /l3(i?i ' ^ + W"^ b/6(i2t ' 

m t m b X bgRL m t m b X t g RL 
24MJ ^h2-2h 5 )(R t ,R b ) __^/ 14 ( jRt)J R 6 ) 

Mgm t m fe X t X b c/ RL m 2 X t g LR 

12Mf /9(jR " ^ + 24M* (/l4 ~ 3 ^^> ^ 

i^ /l5( ^' ^ - lff (m ' /2 + M ^ +/s + m ' /4)( ^' Rb) 



Rh 



M g m 2 X 2 g LR M g m 2 X 2 g LR 

UM$ MRh Rb) 12M| h{Ru Rb) 

M g (M' 2 - M' 2 )g LR M g (M 2 -M 2 R )g LR 
IM, h[Ru Rb) + UMi h(Ru Rb 



(C.15) 
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C.2 Near-zero stop and sbottom mixing 

The results for K t) H t) K b and H b up to order M% w j M 2 SV SY are: 



K t = 



Pis 

3tt 



2M § X t 
' Ml - Mi 

tl *2 



M~ \ i M~ 

In^r + ^(i^tj - gi(R t2 ) - - in— | + 0i (i* 



ti 



m? 

01 



tl, 



- gi(R b i) + 92(Rti) - 92(Rbi) + 



m 2 t Xi 



(Ml -Mir 



1 Mi 

-2 ln W + 9li - Rt ^ ' 9l( - Rt ^ 
ti 



+ g 2 (Rt 2 ) - g 2 (R tl )) 



m\Xl 



(Mi - Mi ) 2 

x hi hi' 



1 M? 

2 + ~~ ^(-^i) + # 2 (#b 2 ) 
61 



mr 



^(i^)) - 6^(2/2 - h)(R bl ) - y^^h - h){Rt > ] 

61 t 2 



^ (2/ a -/,)(*,) + ^"'^ 



if* = 



+ 



12M^ 

»2 

2a, 



Mi - Mi 

01 02 



3n m t cot (3 
M~ g m t X t 



M g g RL 



2Mi 
ti 



f 1 (R t2 ,R l 



In, 



ti 

h(Rbi) 

3Mi 

01 

m t 9LL 
QMi 



3Mi 

02 



(2f w - f 13 )(R tl ,R 



bi, 



+ 



Mi - Mi 

tl *2 

M g m b X b 
Mi - Mi 

bi b 2 



gn + 



m t X t 



gnu 



Mi - Mi 

tl t2 

m b X b 

Mi - Mi 

bi b 2 



gRL 



fi(R t2 ,R bl ) fi(R tl ,R bl ) 



2Mi 

t 2 



2Mi 

ti 



gRL 



f 1 (R t2 ,Rb 2 ) fi(R t2 ,Rb 1 , 



2Mi 

*2 



2Mi 

*2 



m t m b X b g LR 

QMi (Mi - Mi ) 

ti v 61 62 



(fi3(Rti, Rbi) — fn{Rtn Rb 2 ) — 2(fio(Rti, Rbi) — fw(Rti, Rb 2 ))) 



m t m b X t g LR / f 13 (R t2 , R b2 ) fi 3 (R tl , R b2 ) 



Mi - Mi 

tl t2 



6Mi 

12 



6Mi 

ti 



1 m bgRR t /p p \ 

6M? ^ * 2 ' b2 > 

t2 



+ 



fi(Rti, Rb 2 ) fi(Rt 2 , R b i) h(Rti,Rbi) fi(Rt2,Rb 2 ) 



M g m t m b X t X b g LR 

( M ?i- M V( M ii- M V . 

m tX t g R L I fi3(Rti, Rbi) fio{Rti,Rbi) fi3{Rt 2 ,R b i) . fio(Rt 2 ,Rbi 



2Mi 

ti 



2Mi 

t 2 



2Mi 

ti 



2Mi 

12 



Mi - Mi 

tl *2 



6Mi 
ti 



3Mf 
ti 



6Mi 

t 2 



+ 



3M| 

t 2 



m 2 b X b g RL 



QMi (Mi - Mi ) 

52 v 61 62 ' 



(fi3(Rt 2 , R b2 ) - fi3(Rt 2 , RbJ) 



|(ffK72 + m| + / 3 + 7^/4)^2, Rbi) 

i 2 



K b = 



tts 

3tt 



2M g X b 



M 



M| - Mi 

bi b 2 



In 



b 2 



Mi 

bi 



+ gi(Rbi) - gi(Rb 2 ) 



1 Ml 
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- 9i(Rh) + 92(Rb 1 ) - 92(Rt 1 
+ g 2 (R t2 ) - g 2 (R tl )) + 



1 M~ 

{Ml-MlY \ -2 ln Mf + 9l ^ ) ~ 9l{Rt ^ 



mjxl 



(Ml - Ml) 2 \ 2 Ml 

v h b 2 ' \ h 



- ln-£ + gi (R b2 ) - gi (R bl ) + g 2 (R b2 ) 



m 



- 92(Rb 1 )) - ^ri 2 ^ - h)( R ti) - l2M? 

ti 62 



m t (t?f fVR \ 1 M- g m 2 X t 
12M? [2h ~ h)[Kt > } + M~ - M~ 

t 2 ti t'2 



(2f 2 - f 3 )(R b2 ) 

f2(Rh) f2(Rt 2 ] 



3M? 

tl 



3Ml 



H h 



2a, 



37r mfetan (3 



Mgg LR . rn t g RR 

-JlK*Hi, K b 2 ) + O ,,o (^Jio - J 13 ){^t 2 ,rCb 2j 



2M? 

ti 



6M? 

t 2 



M g m t X t ( mtX t \ ( h{R t2 , R b2 ) fi(R tl ,R b2 ) 

+ — 7To~ I fi'RR — TTo TTT9LR 



+ 



9ll + 



M? - Ml 

ti t 2 

m b X b 
Ml - Ml 

bi b 2 



2M? 

12 



2M? 

ti 



9lr 



fi(R tl ,Rb 2 ) fi(R tl ,Rb 1 ] 



2Ml 

ti 



2Ml 

ti 



M? - Ml 

ti t 2 

M g m b X b 
Ml - Ml 

bi b 2 

m t m b X b g RL 
6M? (Ml - Ml ) 

t 2 v bi b 2 ' 

m t m b X t g RL ( fi 3 (R t2 ,R bl ) fi 3 (R tl ,R bl ) \ m b gLL , n , 

+ ^,^0 J^K^h, tibi) 



(fi3(Rt 2 , Rbi) — 2fio(Rt 2 , Rbi) — fn{Rt 2 -,Rb 2 ) + 2fio(Rt 2 , Rb 2 )) 



Ml - Ml 

ti t 2 



6Ml 

t 2 



6Ml 

ti 



6Ml 

ti 



fi(Rti,Rbi) fi(Rti,Rb 2 ) fi(Rt 2 ,Rb 2 ) , fi(Rt 2 ,Rb! 



M g m t m b X t X b g RL 
(Mf -M?)(M? - Ml) 

v ti t 2 ' x hi 62 

m t^t9LR ( fn(Rti, Rb 2 ) fw(Rti,Rb 2 ) fis(Rt 2 ,Rb 2 ) . fw(Rt 2 ,R, 



2Ml 

tl 



2Ml 

tl 



2Ml 

t 2 



+ 



2Ml 

t 2 



6M? 
ti 



Ml - Ml 

ti t 2 

m 2 b X b g LR N 
6Ml(Ml - Ml) 1/131 Ut " Ub *< 

ti ^ bi h "' 



3Ml 
ti 



QMl 

t 2 



+ 



b 2 , 



3M? 

t 2 



fn(Rti, RbJ) 



b 2 > 



^f(rn 2 t f 2 + M 2 H+ h + mlU)(R h ,R b2 ) 



(C.16) 



C.3 Maximal sbottom mixing and near-zero stop mixing 

The results for K t) H t) K b and H b up to order M% w jM 2 slJSY are: 



3tt 



2M- g X t 

'Ml - Ml 

ti t 2 



Ml \ 1 Mf 

ln W + 9l(Rtl) ~ 9l(Rt2) ~ 2 ln W + 9l(Rtl) 



ti 



9i(Rb) + 92(Rti) ~ 92{Rb) 



1 Mf 
, --l n *2. 

(M|-M£) 2 \ 2 Ml 



9i(Rt 2 ) -9i{Rt 



t 2 > 
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+ g 2 (R t2 ) - g 2 (R tl )) 



2 2 
T77 771 



m' 2 - if' 2 

- ^i(4/2 - 2/s + h){R b ) + * * (3/i - f2)(R b ) - 



m 2 Xf 
24M 4 



H, 



2a, 1 



12MJ 



6M| 



(2/3-/4)(i?6) 



+ 



37r m t cot /3 



' RL -fi(Rt 2 ,Rb) + -ft^t(2/io ~~ f 13) (Rh,Rb) 



2M? 

t-2 



M~ — M~ 

tl t-2 



9ll + 



m t X t 



M~ — M~ 
ti 12 



9rl 



6M? 

ti 

f 1 (R t2 ,Rb) fi(Rti,Rb 



2M? 

*2 



2M? 

ti 



+ 



m b g RR 
6M? 

t2 



fn{Rt2-,Rb) + 



Mg9RR 

QMf 

*2 



fn b X b fQ(R t2 , R b ) 



m t m b X b g LR 

24M? 
ti 



(3/12 Rb) — 2/15(^1, Rb)) 



m t m b X t g LR / fi 3 (R t2 ,R b ) f 13 (R tl ,R b 



M? - M~ 

tl *2 



6M? 

12 



6M? 

ti 



+ 



M~ g m t m b X t X b g LR I f 6 (R tl , R b ) f§(R t2 ,R b 



{Ml- Ml) 



<6Mf 

ti 



QMf 

*2 



m\X t g RL ( f 13 (R tl ,R b ) f 10 (R tl ,R b ) f l3 (R t2 ,R b ) f m (R t2 ,R l 



M? — M~ 

tl (2 



6M? 

ti 



3M? 
ti 



6M? 

12 



+ 



3Mf 

12 



!! f^/i5(^2, R b ) ~ ^K 2 / 2 + Ml + f 3 + mlf 4 )(R t2 , R b ) 

h t2 

M~ g m\Xlg RL M g {M' 2 - M' 2 )g RL 

-}s{Kt 2 ,K b ) H 77T774 fa{-Kt 2 ,ttb) 



\2Mf 

t2 



12M? 

*2 



3vr 



X h 1 Af? 

^A(it: 6 ) + - - 9l (R tl ) + 9l (R b ) - g 2 (R tl ) + g 2 (R b ) 

s A M s 



™ 2 X 2 



(Ml -Ml) 2 



1 Mf 



12M? 

*2 



(2/ 2 -/ 3 )(^ 2 )- 



mi 



12M 2 S 



(2f 2 -f 3 )(R b ) 



m 2 X 2 M' 2 - M' 2 

+ it (4/2 - 2A+/4)W) -^W L(3/l - /2)W) 



+ 



H b = 



M g m\X t f 2 (R tl ) f 2 (R, 



m 



M? -M? 

ti t 2 



3M? 

ti 



3M? 

t2 



6M? 
ti 



U2f 2 - f 3 )(R tl ) 



2a a 



1 



37r mfetan (3 



9 ^ 2 h{RtnRb) + J^ 2 R (2/io - fi 3 )(R t2 ,R b ) 



2Mf 

ti 



6M-f 

12 
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M g m t X t / mtX t \ / h(R t2 ,R b ) _ fi(R tl ,R b ) 

+ M?-M?\ 9rR M?-M? 9LR \ 2M? 2M 2 

tl *2 \ *1 *2 / \ *2 *1 

+ f^r~ fi3(Rti,Rb) + rn b X b f 6 (R tl , R b ) 

ti ti 

* 24Af 4 ~ ( 3 ^ 12 ^* 2 ' ^ ~ 2 fi5(Rt 2 ,Rb)) 

_ m t m b X t g RL ( f 13 (R t2 ,R b ) _ f 13 (R tl ,R b ) 
M? - M? \ QM? QM? 

M- g m t m b X t X b g RL ( f 6 (R tl ,R b ) _ f 6 (R t2 ,R b ) 
(Ml -Ml) { QMl 6Mf 2 

m t^t9LR ( fi3(Rh,Rb) fio(Rh,Rb) fn(Rt 2 , Rb) ^_ fio(Rt 2 > Rb) 



M~ — M~ \ 6M? 3M? QM? 3M? 

ti t'2 \ ti ti t'2 t'2 

- ^^MR^Rb) - |g?K/ 2 + Ml + / 3 + mlf 4 )(Rti, R b ) 

tl tl 

Mgrr^QgLR . M g (M' L 2 - M' 2 )g LR /" 
12A0 h{ K ti,rC b ) V2M 1 M-Kti, -K&J 

C.4 Near-zero sbottom mixing and maximal stop mixing 

The results for K t , H t , K b and H b up to order M^ w /Mg USY are: 



M-X 1 M ^ 

-^fi(Rt) + 2 - <?i(^>i) + <?iM - ftiCi^) + <fe(i2t) 

m 2 X 2 / 1 

' ' " ' + 9l (R b2 ) - 9l (R bl ) + g 2 (R b2 ) - <? 2 (i? 6l ) 



(C.17) 



(M, 2 - M- 2 ) 2 \ 2 M- 2 

Ol 02 \ 61 

m 2 m 2 X 2 M 2 — M 2 



2 



M g m\X b ( f 2 {Rbi) hiRbS] ™l m b 
+ Ml -Ml [^M 2 3MT i " 6Mf ( " K J " l2Mp /2 " /sK ^ a) 

01 02 \ 01 02 / 01 62 



2a,, 1 



37T m t cot /3 



fi(R t ,R bl ) + 2 (2/10 — f 13) (Rt, R b i) 



2M| ^^- t '" ^ 1 6M| 



+ m t x tfs{Ru Rbi) + fi3(Rt, Rb 2 ) 

M g m b X b ( m b X b \ f h{R t) R b2 ) fi (Rt, Rbi) 

Ml - Ml \ 9RR Ml - Ml 9RL j V 2M 2 2M 2 

61 62 \ 61 62 / v * * 

m t m b X b g LR 



6M 2 (M 2 -M 2 ) 



(fniRt, Rbi) — 2fi (R t , R bl ) — fniRt, R b2 ) + 2/10 Rbi)) 



30 



m t m b X t g LR M g m t m b X t X b g LR 

71 Ju{ K t, ftb 2 ) + ,. i i , : > 77TT IM-ttt, ttbi) - h{^H, tib 2 )) 



■2\.\I; 

m 2 t X t g R L 
24MJ 



6Ml(Ml-M? 2 ) 



(3/n (Rt, -RbJ — fu(Rt, RbJ) 



m 2 b X b g RL 



QMI(MI-M? 2 ) 



{fn{Rt, Rb 2 ) — fi3{Rt,Rbi)) 



mjh ^Ml + h + miU){R t ,R bl ) - ^r^^fr^RbJ 



M g !hn. u , rl f , 

2AMi 



M g {M 2 -M 2 R )g RL f(jD 
12M| h(Rt,R bl ) 



12Mf 



"M- 2 - Ml 

b\ b'2 



M 



1. M( 



- gi (Rt) + g 2 (R bl ) - g 2 {R t ) + 



m 2 X 2 



(M- 2 - M, 2 ) 2 I 2 M| 



1 Ml 

o /n T7f" + 9i(Rb 2 ) - 9\{Rb 



+ 92(R b2 ) - 92(R bl )) - Y^|r(2/ 2 - fs)(R b2 ) 



b 2 

m 2 m 2 X 2 M 2 — M 2 



H h 



M- g m 2 X t 
2a, 



(2f 3 -f 4 )(R t ) 



+ 



+ 



3n m fe tan (3 
M g g RR 



' 2M 2 ^^ Ru Rb ^ + QM 2 ^ 10 ~~ •^ 13 ^- R *' ^° 2 ) 



6M| 



™>tX t f 5 (Rt, Rb 2 ) + Tf^JT' hsiRt, Rbj) 



6M 2 



M g m b X b 
Ml - Ml 

b\ b 2 



9LL 



m b X b 



Ml - Ml 

bi b 2 



9LR 



( fi(Rt, Rb 2 ) _ fi(Rt, RbJ 
V 2M 2 " 



2M| 



6Mf(M- 2 -M- 2 ) ^ 13 ^*' ~ 2 / 1 °(- R *> ^i) _ /i3(-Rt, ^> 2 ) + 2/i (-Rt, i?6 2 )) 
-fu{Ru RbJ + 



m t m b X t g RL 



M g m t m b X t X b g RL 



24MJ 

m 2 X t g L R 
24Mj 



6MKMJ-M 2 ) 



ihiRuRbt) — h(Rt,Rb 2 )) 



(3fn(R t , Rb 2 ) — fu(Rt, Rb 2 )) 



m 2 b X b g LR 



6M 2 (M 2 -M£) 



(fi3(Rt, Rb 2 ) — fi3(Rt, Rb 1 )) 



M g m\Xtg LR 



M "- lLR <m\h + M 2 H+ f 3 + m 2 b f 4 )(R t ,R b2 ) - fr(Rt, Rb 2 ) 



24M% 
M g (M 2 -M 2 R )g LR 
+ 12M| HRt,R b2 ) 



(C.18) 
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